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ABSTRACT 


The various wave theories, theoretical breaking criteria 
and derived breaking criteria are reviewed for shallow water 
waves. To account for the non-linear hydrodynamics present. 
in a shallow water wave breaking on a beach with a sloping 
bottom, the perturbation technique of Iwagaki and Sakai is 
used to derive a second order expression for the horizontal 
water particle velocity for long waves. The kinematic 
breaking criterion is applied to the derived c(2) and ul2) 
values to establish breaking. The results indicate that 
the ratios of 7,/L 5 and h,/H, provide reliable breaking 
criteria. Each of the parameters is dependent only upon 
beach slope and H,/L,. Theoretically derived values for 
h,/H, compare favorably with field measurements and offer 
Significant improvement over previous theory. Predicted 
breaking depths are less than those present in experimental 


data, suggesting extension to higher orders may be warranted. 
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I. INTRODUCTION 


Attempts at deriving breaking criteria have been made 
since Stokes (1847) presented his classical development. 
Derivations of the many available water wave theories all 
involve the solution of Euler's equations of motion coupled 
with the continuity equation for incompressible, inviscid, 
meeotational flow subject to certain boundary conditions. 
Breaking, or near breaking waves have very steep profiles 
in which the wave height is large compared to the relevant 
length scale implying the hydrodynamics are highly non- 
imdear. At the onset of breaking strong vorticity is 
introduced at the surface near the crest and the assumption 
that the motion is irrotational is no longer valid. The 
strong abe ac and induced vorticity make the 
analysis of breaking waves mar haat teats Veryraqitricul ce. 
This thesis is concerned with finding an incipient breaking 
criterion for waves shoaling on a beach before vorticity 
is induced but including nonlinear effects. The discussion 
will be limited to shallow water wave theory. 

Solution of the equations of motion has required the 
application of physical assumptions associated with various 


wave characteristics. Thus, each formulated theory is 
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limited in its range of applicability to regions where 

its underlying assumptions are valid. Shallow water wave 
theory may be classified according to the bottom being 
horizontal or sloping. This is an important restriction 
because field and laboratory measurements of breaking waves 
suggest that the bottom slope is an important parameter in 
the classification of breaking waves. 

The simplest form of solution to the wave problem is 
to linearize the equations of motion assuming a horizontal 
bottom. Peregrine (1972) has shown that the linearized 
equations apply to regions where the ratio H/L and H/h are 
both much less than one, where H is wave height, L is wave 
length, and h is water depth. Consequently, this solution 
is restricted to waves of infintesimal height and not 
applicable to steep breaking waves. 

Stokes (1847) was the first investigator to present a 
higher order solution applicable to finite amplitude waves, 
though limited to a horizontal bottom. In his development, 
he transforms the basic equations to an equivalent set by 
using a velocity potential, ¢ . The solution is obtained 
by expanding the velocity potential using a perturbation 
scheme which employs H/L as a perturbation parameter. To 
the lowest order, Stokes’ method results in linear theory. 


The accuracy of Stokes’ solution at a particular order 
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decreases as the ratio h/L decreases. Dean and Eagleson 
(1966) attribute this inaccuracy to increasing bottom 
influences and a decrease in the importance of vertical 
particle acceleration. De (1955) concluded that Stokian 
theory should be discarded for values of h/L of 0.125 
and less. Dean (1968) expanded a stream function using 
a numerical perturbation scheme and was able to raise a 
"Stokes'' type wave to any desired order. In this manner 
he was able to solve for incipient breaking. 

When the relative depth is very small, as in very shal- 
low water, the vertical acceleration can be neglected and 
the fluid path curvature is small. Hence, the pressure 
is assumed to be hydrostatic as the vertical component 
of motion does not influence the pressure distribution. 
The resulting equations are referred to as the "long wave 


equations." 


A sloping bottom and finite amplitude are 
allowed by the long wave equations. 

Freidrichs (1948) was able to derive the long wave 
equations by a rigorous mathematical approach. Utilizing 
quantities h and L which represent typical depth and length 
scales, the Airy equations were non-dimensionalized. This 
procedure resulted in a large stretching of the horizontal 


coordinate relative to the depth coordinate. A perturbation 


analysis was then applied with the perturbation parameter 
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g = n*/L2, He found that to the lowest order the pressure 
was indeed hydrostatic and that the long wave equations 
resulted. 

A third length scale for shallow water theory, which 


utilizes the wave amplitude a, is the Ursell parameter 
3 


t+ | 
alr 


Ursell (1953) showed that long wave theory belongs to the 


regime 


>So 


tt |m 
aie 


Boussinesq (1872) assumed that the pressure was no 
longer hydrostatic, which allows inclusion of vertical water 
particle velocities, but results in a limitation on the 


wave height. The Boussinesq equations apply when 


3 


| 
oa la 
t= 


implying the waves are not as high and the water is 
relatively deeper than for long wave theory. 

Korteweg and de Vries (1895) simplified the Boussinesq 
equation by considering waves which travel only in one 
direction over a horizontal bottom. Extending Boussinesq's 
equation in this manner they produced a wave theory they 
termed "cnoidal."' The limiting case of cnoidal theory is 


the solitary wave. The cnoidal/solitary theory has received 
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considerable attention by researchers in recent years. 
Keller (1948) extended the perturbation analysis employed 
by Freidrichs to the first approximation. He found that 
to the first order his results were those of Korteweg and 
de Vries. Laitone (1960) continued the process further, 
obtaining second order approximations to cnoidal/solitary 
waves by solving Freidrichs’ method to the fourth order. 
The assumption of waves traveling only in one direction 
precludes a reflected wave and imposes the important 
restriction of a horizontal or nearly horizontal bottom. 

The long wave equations 2622 used in this study because 
it is felt that properly including the sloping bottom is 
the most important next step in seeking a breaking wave 
criterion. The possible importance of vertical accelera- 
tions in the wave breaking process are recognized, but 
are assumed negligible in order to obtain mathematical 
tractability. A second order solution of the long waves 
is sought and a breaking criterion derived based on a 


kinematic instability condition. 
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II. THEORETICAL BREAKING CRITERIA 


It is desired to formulate breaking criteria which can 
be expressed in terms easily observable and measurable. 

The several breaking criteria which have been developed 
may generally be broken into two categories. First, there 
are those which are derived from waves of steady form. The 
waves considered are assumed to be in shallow water of 
constant depth. Both Stokian and cnoidal/solitary wave 
theories have been employed in these investigations. The 
second group of derived criteria consists of those con- 
cerned with waves which deform as they shoal. Research in 
this category of waves has been confined to long wave 
theory. The derived criteria for horizontal and sloping 
bottoms will be investigated separately. 

In order to determine wave parameters which can be used 
to predict the breaking of waves, it is necessary to first 
formulate some type of limiting value. The physically 
significant breaking criteria are the kinematic, dynamic 


and geometric criteria. 
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A. KINEMATIC BREAKING CRITERION 

Originally formulated by Rankine (1864), the kinematic 
breaking criterion states that the limiting value of the 
water particle velocity at the wave crest, u,, is the wave 
fase velocity, ¢, u.=c. Physically, this is a logical 
limitation, for if the particle velocity is allowed to 


exceed the wave phase speed, the particle would separate 


itself from the wave form. 


B. DYNAMIC BREAKING CRITERION 
The dynamic breaking criterion is stated by the vertical 


momentum flux equation at the surface (z=7”) 


Diwan OD - 
Dt p dz S: 


The criterion can be formulated in several ways. Assuming 


the pressure is a constant at the surface, the maximum 


acceleration is 


ae - p constant 


If the water particle acceleration exceeds the gravitational 
acceleration, the water particles leave the surface vertically. 

A second statement of the dynamic criterion concerns 
the vertical pressure gradient. Laitone (1963) proposed 
that the limiting vertical pressure gradient beneath the 


wave crest is zero. In his study of cnoidal and solitary 
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waves he determined a value for H/h beyond which the 
pressure gradient reverses its sign and the pressure begins 

to decrease with depth. He concludes that this is physically 
impossible and thus accepts zero as the limiting pressure 


gradient value. 


C. GEOMETRIC BREAKING CRITERION 

As a wave progresses into shallow water, the surface 
slope steepens. The geometric breaking criterion places 
a limiting value of infinity (vertical face) on the slope 
of the water surface. Beyond this value, the wave becomes 
unstable and the water particles fall forward ahead of the 
wave, Stoker (1957) shows that a vertical slope can be 
obtained. An insight into the concept can be obtained by 
considering the speed of the shallow water wave disturbance 
to be given by c = [a(n th)]2. Since the crest of the 
wave has a greater depth of water beneath it than the 
trough in front of it, it tends to "catch-up'' with the 
trough. Hence, the forward face continues to steepen as 
the wave shoals until it is vertical. 

A review of the application of the three breaking 
criterion to the separate categories of waves of permanent 


form and shoaling waves follows. 
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III. DERIVED BREAKING CRITERIA 
A. BREAKING WAVES OF PERMANENT FORM (SHALLOW WATER, 

CONSTANT DEPTH) 

The theoretical kinematic and dynamic breaking criteria 
have been applied to waves of permanent form. The former 
has resulted in a number of derived criteria while the 
latter has only been applied to cnoidal/solitary theory to 
produce single breaking parameters. 

l. Kinematic Breaking Criterion 

a. Crest Angle 

By applying the kinematic breaking criterion 
to his formulated wave theory, Stokes (1880) produced the 
first derived breaking criteria. He showed that when the 
enclosed crest angle, Figure 2, reached 4 x radians, 
120°) , breaking would occur. To arrive at this conclusion 
Stokes was forced to make two assumptions: 1) that the 
crest would be formed by two intersecting straight lines 
tangent to the real water surface curvature, and 2) that 
the velocity potential, transformed to polar coordinates, 
could be approximated by 

b(r,¢@) = Br sin (NB) , 
where B and N are to be evaluated and r and @ are polar 


coordinates. 
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ENCLOSED CREST ANGLE FOR KINEMATICALLY LIMITED WAVE 


Figure 2 
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Several investigators have verified Stokes' 
criteria (Gaughan, Komar and Nath, 1973), however the range 
of depths over which Stokian theory can be applied places 
limitations on this criteria. Therefore, a contribution made 
by Chappelear (1959) in which he was able to develop a method 
of satisfying the kinematic breaking criterion for all depths 
is of particular importance. Chappelear also verified 
Stokes' value for the limiting crest angle. 

b. Wave Steepness 

Another breaking criteria which has gained 
acceptance is that of wave steepness in which it is stated 
that the wave height is limited to one-seventh of the wave 
length. Several researchers, including Michell (1893), 
Havelock (1918), Davies (1952), Yamada (1957) and Chappelear 
(1959), have produced results close to this figure. Their 
values for maximum wave steepness vary from 0.1412 to 0.1443. 
All of their derivations closely followed classical lines. 
Dean (1968) used a numerical stream function approach to 
examine breaking wave criteria. He found a solution to 
the full nonlinear wave formulation that is exact except 
for the dynamic free surface boundary conditions. The 
solution is found by fitting the dynamic free surface 
boundary conditions numerically using an iterative scheme. 


Taking the limiting horizontal particle velocity as 98.5% 
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of the wave phase velocity, the resulting maximum wave 
steepness value was 0.1723, which differs from those 
previously given by about 20%. 

c. Wave Height to Depth Ratio 

The kinematic breaking criterion has been used 

by many investigators to obtain limiting values for the 
wave height to depth ratio, H/h. The results, however, 
have not been at all consistent and vary widely with 
different wave theories. Chappelear (1959) analyzed 
Stokian waves and arrived at a value of 0.87. For solitary 
waves, results obtained by Boussinesq (1871), Rayleigh 
(1876) , McCowan (1894), Gwyther (1900), Packhan (1952), 
Davies (1952), Yamada (1957), Lenau (1966) and Yamada, 
Kimura and Okabe (1968), vary from 0.73 to 1.03. Dean (1968) 
found a value of 1.0 using his numerical stream function 
approach. Gaughan, Komar and Nath (1973) express the 
opinion that these differences probably arise due to 
approximate fits of the complex velocity potential to the 
free surface boundary conditions. As will be discussed 
in a later section, Laitone (1963) found values of H/h = 
0.73 and 0.81 for solitary waves using different theoretical 


breaking criterion. 





2. Dynamic Breaking Criterion 
a. Water Particle Acceleration 
Various Limiting values of water particle 

acceleration have been determined. Kinsman (1965) gives the 
gravitational acceleration, g, as the limiting value for a 
crest angle of 120°. Gaughan, Komar and Nath (1973), how- 
ever, show that by using Stokes' wave crest equations and a 
crest angle of 120s8 the limiting value should be &, 
b. Vertical Particle Velocity 

Laitone (1963) examined the vertical water 


particle velocity for cnoidal and solitary waves. Using 


a third order velocity equation for cnoidal waves, he found 


8° 
90 +2 


erty of the wave, is restricted to 5<@41, produced 


that values of H/h greater than , where ©, a prop- 





physically impossible velocities. Hence this ratio was 
established as the limiting value. As the value of 6 
approaches the limit of 1.0 the cnoidal wave approaches a 
solitary wave form. Therefore, the limiting value for a 
solitary wave to the third order is H/h = = 0.7272. 
Laitone carried his solution for the solitary wave to the 
next higher order and found a value of H/h = 3-15 


0.7321. These two results compare favorably. 
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c. Vertical Pressure Gradient 
Laitone (1963) developed a different limiting 

value for H/h than that previously discussed when he applied 
the vertical pressure gradient criterion. In this case he 
found that for H/h = (2 @ /3)%, 5<@<1, the pressure 
gradient is zero for cnoidal waves. This expression was 
derived to the third order. If the ratio H/h increased 
from this value, the sign of the gradient would reverse, a 
condition that he concluded could not exist. [In the limit- 
ing case of @ = 1.0, which gives a solitary wave, the 
limiting value of H/h = 0.812 is obtained. This differs 


significantly from his previous limit of 0.7272. 


B. WAVES WHICH DEFORM AS THEY SHOAL (SLOPING BOTTOM) 

As a wave shoals over a sloping bottom, the wave height 
and profile are altered. The theory of long waves has been 
most generally applied to research in this region. The 
derived breaking criteria have been formulated through the 
use of the kinematic and geometric breaking criteria. 

l. Kinematic Breaking Criterion 

The first attempt at applying the kinematic breaking 
criterion to a deforming wave was made by Ayyar (1970). His 
derivation made use of the concept of a wave front. Simply 


stated, a wave front is the position where a discontinuity 
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in the surface slope occurs. The slope of the surface 
will be zero in front of the wave front and negative behind 
it. Ayyar's approach was to obtain the slope at the wave 
front, integrate to find the free surface # , and then 
apply the kinematic breaking criterion. Assuming the 
geometry of the plunging breaker, he then derived the 
value of Yp/hy, =e enmwienre sy pea emcight of breaker crest 
above the bottom and h, = depth at breaking point below 
the still water level. 

Several problems exist in Ayyar's derived criterion. 
First, the derivation is based on the geometry of a plunging 
breaker and thus excludes the other categories of breaking 
waves. Additionally, his formulation assumes that breaking 
will occur at the wave front. Gaughan, Komar and Nath 
(1973) point out that this may not be a valid assumption. 

2. Geometric Breaking Criterion 

Use of the geometric criterion has been made by 
several researchers in developing long wave breaking 
criteria. Stoker (1957) showed that long waves could 
obtain a surface slope of infinity. He extended his work 
and used a numerical methods technique to arrive at a 
solution to the problem. The method, however, is not 
satisfying in that it requires a number of approximations 


to be made and requires recalculation as the initial 


26 





conditions are altered. Further discussion of this pro- 
cedure will be presented in a subsequent section. Burger 
(1967) and Greenspan (1958) used the wave front concept 
and the vertical slope criterion to predict the horizontal 
distance traveled from the wave front at time t=0 to the 


point of breaking. The result was 


aun (2 2s.) 
S+M 


puieres = slope of the beach, h(x) = hy-Mx, h), = initial 
water depth, and S = initial surface slope at the wave 
front. Gaughan, Komar and Nath (1973) discuss several 
limitations to this result. As was the case in Ayyar's 
work, the breaking is assumed to occur at the wave front. 
The surface slope behind the front is not examined. Some 
other point, such as the wave crest, may become vertical 
prior to this condition occurring at the wave front. Also, 
the use of horizontal distance to breaking is not a useful 
criteria. It is not easily measurable, having a somewhat 
arbitrary origin. Prediction of the wave height at breaking 
is a much more useful parameter. 

Another approach involving the vertical surface 
slope was taken by Mei (1966). The basis for his work was 
originally proposed by Carrier and Greenspan (1959). The 


technique is to produce a set of characteristic equations 
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from the long wave equations and then to make an additional 
transformation through the use of a final pair of inde- 
pendent variables, o and ». The equations are reduced 
to a single linear equation which involves a velocity 
potential, o( ) 3 

(so), - T$,,=0 
Mei solves this equation to the first order. He then 
follows a procedure used by Carrier and Greenspan whereby 
the Jacobian J = 0d(x,t)/ a(¢,%) is investigated. This 
Jacobian will vanish at points for which the surface slope 
is infinite. Mei was able to obtain an expression for 
hi, /H, which was dependent upon H,, L, and the bottom slope. 


The subscript , denotes deep water conditions. Unfortunately, 


O 
Mei found that his theory compared poorly with experimental 
data. Predicted breaking depths were too large, dependence 
on the beach slope was too great and the wave profiles were 
too sinusoidal. Mei suggests that these discrepancies 


could possibly be eliminated by extending the solution to 


a higher order. 
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IV. REVIEW OF LONG WAVE SOLUTIONS 


Essentially two approaches exist which can be used to 
solve the nonlinear long wave equations. The first pro- 
cedure, initially formulated by Stoker (1947), makes use 
of a solution technique known as the method of character- 
istics. A final solution utilizing this method can be 
made either through numerical calculation or by an analyt- 
ical approach. Iwagaki and Sakai (1972) propose a second 
Peiiit ion procedure which involves a perturbation expansion. 


An evaluation of each of these techniques follows. 


A. METHOD OF CHARACTERISTICS 

The application of the method of characteristics to the 
long wave equations was explained by Stoker (1958). Neces- 
sary to this development is the acceptance of the wave 
phase speed relation c = (g( th)]%. The validity of this 
equation is discussed in a subsequent section. From this 


expression it is seen that 


Cx = (gn, + gh) /2¢ oe 


and 


C42 97, /2c. (2) 
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The long wave equations are 

ee EON. x (3) 
and 

Lu(y+h)l,= "ee (4) 
Substitution of (1) and (2) into these equations yields 

Up tuuy t2ec,—ghy=O (3a) 


and 


ac ,+2uc, + cu, =O (4a) 


Relations (5) and (6) result from the respective addition 


and subtraction of (3a) and (4a), 
{ a/At + (ure) ox} Cu+2c-gh,t)= O (5) 


{fat + (ure) /ax}-(u-2e-gh,t)= @ ne) 


The interpretation of (5) and (6) is essential to the 
development. (5) implies that the function (ut2c-gh,t) 
remains constant for a particle moving with a velocity of 
ute. A similar evaluation of (6) can be made. In other 
words, two characteristic curves, C, and Cy, are defined 


such that 
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Ci: dx/dt= ure 

and : (7) 
Co: dx/dt=u-c 

where | 


u+2c-gh t= K = A CONSTANT ALONG on 
and (8) 
UuU- 2C- 9 hb, t= K,> A CONSTANT ALONG Cs 


The system of equations given by (7) and (8) is equivalent 
to that defined by (3a) and (4a), hence, a solution of either 
set provides a solution for the other. 

lt now becomes necessary to determine the character- 
istics so that a solution may be obtained. Stoker's 
technique was to make use of a method of successive approxi- 
mations. He assumed that the initial values of u and 77 
could be observed. Use ofc = Ce(77 +h)? could then pro- 
vide the phase speed at time t=0. These initial conditions 


are assigned values such that for t=0 


u(x,0) = i (x) 
c(x,0) = €(X) 


The task is to approximate u and c for small increments of 


(9) 


time. Figure 3 clarifies the discussion. 
A series of points along the x-axis, which are separated 
by a small distance db,» are considered. Since the values 


of u and c are known for each of these points from (9), the 


Slope of the characteristics C, and C, at each point can 
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INTEGRATION BY FINITE DIFFERENCES . 


Figure 3 
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also be obtained through the use of (7). These slopes are 
used to construct straight line segments from the points 
along the x-axis. Location of points 5, 6 and 7 is deter- 
mined by the intersection of these line segments. A source 
of error is inherent in the use of straight line segments 
to approximate the characteristic curves. This error is 
restricted to a minimum by using Smee evenery small values 


of § For this case, the tangents to the curves provided 


Be 
by the slopes give good approximations to small segments 
of the curves. Equations (8) and (9) give the character- 


istics issuing from points on the x-axis as 


along C, : u+2e-gh,t=ur2c 
(10) 

along C5 : cde Ss ue 
The values of x and t can be obtained for the points 5, 6 
and 7 (this could be accomplished graphically, for instance). 
Equation (10) can then be used to determine the values of 
u and c at these points. The procedure can be continued to 
obtain values for u and c at the subsequent points 8, 9 and 
10. In this manner, a net of points at which values of u 
and c are approximated could be constructed which covered 
an entire field of concern. Values at intermediate points 


could be found by interpolation. Stoker states that as 
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6 0, the process will converge to an unique solution 
of (7) and (8). 

Numerous objections have been raised against this 
particular numerical approach. Ayyar (1970) categorizes 
these into three areas. First, the solution is not explicit 
and requires re-calculation for a change in initial condi- 
tions. This alone makes the procedure difficult to use 
effectively. A second criticism is that a spilling type 
breaker is always predicted; thus the elimination of the 
several other breaker types places a severe Peeenerion on 
the method. Finally, LeMéhauté (1968) has given evidence 
that the technique produces incorrect predictions of breaker 
points. 

An alternative and more analytical solution to equations 
(7) and (8) has been offered by Carrier and Greenspan (1958) 
and was discussed in Section III.B.2. Two independent 
variables ¢ and A are used to transform the equations. 
Defining a Pieity potential ¢(¢,A) yields the linear 


second order equation 
(66) - 0%... = Or aia) 
to which the authors propose a solution, 


o(¢ d)= Ad. (uc) cos( way) (12) 
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where A = constant to be determined and ¥= a phase lag. 
At this point, Carrier and Greenspan made use of the verti- 
cal surface slope criterion and determined values of a 
Jacobian, J = 0(x,t)/d(¢ ,A), which would exist for the 
specific case of non-breaking waves. Consequently, the 
remainder of eheir study is of little value to our discussion. 
Mei (1966) considered the solution technique of Carrier 


and Greenspan as it applied to the case of breaking waves. 


He selects 


v=-BLJ,(0/2i%) cos(/2i 24/4) 
+No(9/2i™) sin(rf2i®em/4y | 23) 


where B =a constant to be determined, ( = Slope of the 
beach, J, = zero-order Bessel function, and N, = zero-order 
Neumann (Weber) function as his solution. Using this rela- 
tion, Mei derived first order equations for u and 7{ . The 
coefficient B is determined by matching the solution to 
that for an outer region where a horizontal bottom Airy 
Theory was applied. Mei determined a breaking parameter 
by applying the geometric breaking criterion. 

Comparison with experimental data indicated that Mei's 
results were not accurate. He attributed these deficiencies 
to the fact that the solution was not carried to higher 


orders. Previous studies conducted by Benney (1966) had 
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suggested higher order derivations were required for shallow 
water breaking conditions. Solution to a higher order is 
required in both the inner and outer regions if the coeffi- 
cient B is to be accurately determined. The procedure of 
Tlapa, Mei and Eagleson (1966) provides a perturbation 
expansion for the outer region. Review of this method 
shows that the third order solution must be considered to 
uniquely determine the second order coefficients. Assuming 
that a successful solution to the outer region expansion 
could be found, or that B can be found independent of the 
off-shore regime, difficulties still remain in the near 
shore area. This is the location of breaking and thus of 
concern. The solution would be vastly simplified if a 
method could be found in which only the velocity potential 
used by Mei need be perturbated. Unfortunately, the single 
equation for the velocity potential prevents this approach. 
The quantities 6 and A are only variables used to trans- 
form the characteristic equations and hence perturbations 
on them produce no new equations. The only remaining 
alternative is to use the next higher order expression of 
the Airy equations to form the characteristic equations. 


This implies the consideration of cnoidal theory. 
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B. IWAGAKI AND SAKAI PERTURBATION TECHNIQUE 

Iwagaki and Sakai (1972) proposed a second method of 
solution which involved a perturbation expansion of the 
long wave equations. The premise of their study was that 
the asymmetric profile of shallow water waves and the effect 
of the beach slope on wave transformations could be explained 
by taking into account the nonlinearity of the long wave 
equations. To show this, they developed a second order 
solution for 77 from the long wave equations (3) and (4). 

In the derivation, the beach slope is represented by i 
and the depth, h, is such that h=h(x)=i°x. Iwagaki and 


sakai assume that u and 72 may be expressed as 


nex + x= 7 (2), c (aire) 


u=& UD 4 x%y@ 4 - (15) 


where & is a small quantity and the superscripts (1) and 
(2) indicate the order of the term. Substitution of (14) 
and (15) into the long wave equations, and arranging them 


with respect to « and ae yields the four equations, 


ms Oran, =e) (16) 


7, tLe®h],=0 (17) 


® (a) (2 (2) 
uw, tu uxt 97, =O (18) 
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ay [ wn? + U uh] = (19) 


Equations (16) and (17) together give two equations in two 
unknowns which can be solved to obtain first order expression 
for 7? and u of the form 


7 (x, t)=qLeosct: J. (26V45% jt se ) 
~ sing t-N 9 (20V Feet Zz) | (20) 


and 
u(x,t) =a x2 El wact-d., (20 VE ) 
+cos at-N ANZ ZY] (21) 


where Jf ) and NC ) represent Bessel and Neumann (Weber) 
functions respectively and a is a constant related to the 
wave height. Equations (20) and (21) can be substituted 
into (18) and (19), reducing the problem to solution of 
(18) and (19) for second order expressions for u and 7 . 


(2) 


The authors presented a solution for 77 in which they 
approximated the Bessel and Neumann functions in the first 
order solutions with trigonometric functions. They show 


that for values of x greater than ci/ag? these substitutions 


are valid. The derived formula for 702) (x, fe): «is 
7x) = Tito: “cos fet +26(3°) "94 
4 x fans AC; a y- 


+ (yy X cos os(3 a ™%) 
= et x sin eee “11 1. i 


16s 
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Iwagaki and Sakai presented experimental results which 
confirmed that their solution technique produced valid 


results for the wave form. 


C. PHASE SPEED RELATION 

In several of the previous sections mention has been 
made of the phase speed relation, c = Ce(7 th] ®. Although 
this equation is widely used in the literature on wave 
theories and is generally accepted; few discussions have 
been presented which establish its validity. The question 
deserves some attention prior to application of the kinematic 
breaking criterion. 

The first inclination toward the use of the relation 
comes from the linearized long wave equations for water of 
constant depth. Stoker (1958) showed that for these 
conditions u satisfies the relation 

Hixx” OF Wee =O; 
7? can be shown to satisfy the same equation. The implica- 
tion of this expression is that the speed of the wave dis- 
turbance is given by c = Gui 

Another indication that the expression may indeed be 
valid comes from what Stoker terms the gas dynamics analogy. 
Stoker credits the development of this idea to Riabouchinsky 


(1932). Consideration is made of a mass per unit area 
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expressed by 
@ = el7+h) (23) 
where Qe = density of water. Thus 


yes C774 ; Cs 


The force P per unit width is defined as 


P =|" pde : (25) 


By using the hydrostatic pressure relation, P may be 


reduced to 
page (y+h) = %e pr. 2 
Multiplying both sides of equation (3) by e (+h) gives 
e (+h) (u,+uu,) = a e (n+ h) 7, (27) 
which may be re-expressed, using (23) and (26), as 
@(uptuuJ=—pytg @hy- (28) 
Equations (23) and (24) can be used to re-write (4) as 
(Zu), =- Sy (29) 


Equations (26), (28) and (29), when combined, give results 


similar to those of gas dynamics for one dimensional flow, 
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the only difference being the presence of the term ge he 
in (28). For the case of constant depth, this term 
vanishes. In gas dynamics, the sound speed is given by 
Cc -L4451* Applying this with feat and (26), , 
a =[ Fs 1*= (538. \* = Z L 3(n+h)] /2. 

Although these aie seine provide an insight into 
using the phase speed relation, they can hardly be consid- 
ered as a definitive argument. A satisfactory derivation 
can be obtained by returning to the methods of character- 
istics. The review of this technique explained that the 
phase speed was assumed to be given by c = Ce(7 +h)] 2. 

The theory could have been derived equally as easily by 
simply defining a quantity c = [g(7 +h)] @. No physical 
significance need be immediately applied to c. This being 
the case, an identical development can be made resulting 

in the same characteristic equations. The task then becomes 
to discover the physical meaning of c. 

Stoker (1958) presents the following argument. It is 
assumed that the initial values of u and 77 are given for 
a body of water which is in motion. The value of c for 
this instant is given through c = Te( seth) 2. Figure 4 
will aid in the explanation. Consider a disturbance created 
over the segment of the x-axis Qi Qo. How will this effect 


the solution? Each point Q on the x-axis has what is termed 
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PROPAGATION OF DISTURBANCES 


Figure 4 
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as a range of influence. This is the region of the x,t 
plane in which the values of u and c are influenced by 

the initial conditions at Q. This area is defined by the 
characteristics issuing from Q. Consequently, for the 
particular case under consideration, the segment Qi Q5 
will have an influence on u and ec for the shaded region in 
Figure 4. The two curves are given by Cy:d,/d, =ute 

and Co:d,/d, =u -c. u is defined as the horizontal 
velocity of the moving fluid. The speed of the disturbance 
as it moves through the flowing water must therefore be- 
given by c. Thus, the validity of the phase speed relation 
seems apparent. 

One final argument can be formulated by the use of the 
method of characteristics. Although he does not discuss this 
aspect of his study, Greenspan (1958) outlines this proof. 
Again, consider the method of characteristics to be formed 
using a term c defined as c= Ce(77 +h)]?. A wave which is 
progressing into quiescent water is considered. Of concern 
is the forward moving characteristic curve, < =) ey 
which contains the wave front. Since the water immediately 
preceding the wave front is quiescent, the value of u must 
be zero at the front. The characteristic, and consequently 


the wave front, must be progressing with speed c. Hence, 


L 
the phase speed is given by c =[.g(7th)]?. 
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V. DERIVATION OF WAVE-INDUCED VELOCITY AND 
WAVE PHASE SPEED USING IWAGAKTI AND SAKAT 


PERTURBATION TECHNIQUE 
A. DERIVATION OF uC?) 

The second order wave induced velocity u62) is derived 
using the perturbation technique used by Iwagaki and Sakai 
to derive a second order surface profile. All terms used 
are defined in the list of symbols. The initial equations 


employed are (18) and (19), 


UL , a + ioe wf Ay n> a 
¥(e) (4) a ' 25 a 
(te h 7 


These equations result when u and 7 are given as power 
series expansions of a small quantity oa and these expres- 
Sions substituted into the long wave equations. Eliminating 
7h?) from (18) and (19), substituting h = i*x and grouping 


terms of u() yields 


(2) : (2) eae 
Mee 9 tix Uy ting f= 


L G@) @) (1) 
Ae 07 ri stad +u 
gt ae Oe ext (30) 
ae uy - Ke 


Using the first order term for u given by iivagaki and Sakai 


and differentiating with respect to x gives 


es 





L(x)=— % (AY 4 fsmot-d, we 


+cos gt - Ni Lp X 


+ Fy fowot (JL pO] HH AL 200]) 
5 coset (iJp(i]- 355 NL Gi) (31) 
where p (x) ea Se (2)*. 


Similarly, differentiating (31) with respect to x provides 


an expression for c-. 


a es N= at (4 a“; 72 3 sine J Leell 


+cosort: N,L P (x) 
~ 25 § swat (dl p6- ty Tp Ge] 
+ cos ie ( Nip (| os N ls I} 


— * aa X “35 wort- dle] 
’ + ott N L (x) 1 §. (32) 


An equation for 7 can be obtained from (31) by 


differentiating with respect to t, 


us . (x, t)= een (a)*x *§ ent Lp] 
—-SINTt.: N iL p(x) J § 


+ ee wae (J. Ep&))-5% AEC) 


a srt (N Cpls] NLpsll) , 633) 


45 





When le) is differentiated with respect to t we obtain 
Y- os 
= (x)= om (2)" 0 =} coset: “JL o(x) | 
=; sme t-N, [a(x 1 > (34) 


Similarly, differentiating 7) (x,t) yields a set of 


equations, 


(1) ma ai’ ] 
(x,t) = eine coset d as (x) | 
P ee oe OF 


and 


ext) = 40g)" x cos v te) lpo] 
(xx = yt iL p@) 13 
= es oo re (3, Cp Calpe stp6o) 
sive t(NeLp ODI Nulpbll)§, 28 


Equations (31) through (36) can be used to find 
expressions for the individual terms on the right hand side 
of (30). When these expressions are combined, the right 


side of (30) is evaluated as 


cos2et 2 Ge i : ame) | PON, Lp (| 
+ J CpOd] No Pay] ‘y 


4+ 96 oe ( J eG] N lp] 
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-2 4, L(x) J N, Lp (x) J) 
~3obe* (4) «(J EpGi] N00] 
care LPG) IN, fo (x)])§ 


+Sinect to! ax MC JC CO]LL (4) 
ies NoLp(x) JN, PO] 
MOF 2 (Ly Posl-J7 [poo] 
Ca etl =] *Tolxdl)f 


+2q(2) Wacsater Lo) IN {Lol i) | =) EL GQ)N, Lal. (37) 


The asymptotic expansions of the Bessel and Neumann (Weber) 


functions are 


J, Cw) ~ ~ (24) eae ZB g-/, \ 
w)~ a oe ( w- r “os — 


the approximations being valid for values of w such that 
[wf =O) 
Substitution of (38) into (37) gives as a final expres- 


sion for the right side of (30), 
(1) (2) (1)_ &) ies (3) a “) W) WW) 
q (zu, Nx tA Myx~ Tx PPO HET EM 


Cos Zot} = 8 X es P(x) | 
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— 15/4. oa (2 \*, Fey [2 P(x) | 
as ane os| 2 s 


+ sin2ot {seat swl2 963 pen Er 8 cosl2gh)| 


ea x “sw l2(n) | 
2 A 3 
ae. x (39) 


neve POE te (S\*- 17 
where P (x)= 20 (9% 4 





The solution of u(2) is assumed to be of the form 


Plxd)= Ald) cos2et+B(x)sin2et +C (x), (40) 


Use of this expression yields equations for the individual 


terms on the left side of (30), 


Lx t) = =A (x)cos2ot-+ BY x) SiIn2ot +C (x), (41) 


1 (es) =A'(x)cos2et+B (x\sin2ot +c"(x) (42) 


UP (x, = 20 { B(x) Gear A(x) sida { (43) 


) 


a” (x t) = -4¢ f(x) coslat + B(x) sindot § (44) 


Combination of (41) through (44) can be used to provide an 


expression for the left side of (30), 
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(2) 


et 7 
Cos2at {- ta" A(x)~gix M6) ee g ( A(x) § 


+sineot §- Ar "B(x)=q xB “(x)- dq (x) § 
: {2 “gx “(x) - zai (x) - (45) 


g (ix Pied = EG Ut . \= 


Comparison of (39) and (45) allows a determination of 
A(x), B(x) and C(x) to be made. This is accomplished for 
A(x) and B(x) by equating the respective coefficients of 
the sine and cosine terms in the two equations while C(x) 
is found by equating the terms independent of time. This 


procedure gives 


A(x) = a (2) atid] 2 P(x) | 


PAN 


sts - Xx = cos [2 (>|, (46) 


B()= 4 (A) *% cosl2 p(s] 
— We 23. een (2 ‘90, (47) 


and 


C(x) = a. ne (48) 


21 oS 


Substitution of (46), (47) and (48) into the left side 


of (30) yields 
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cos2et | “wee x* cos[2 P(x) | 
-B a (% 12 xX ssl 2 P(x) ] 
— of. eee X-F cos LA P(x)] $ 


Te sina tj ee Ks inl2 9x) 
as 'S/, a9 (%)% /2 ee f2 P(x) | 


+ &fs og 3 sin x 
= [299 1 § 








Saas, 
4 .K* . (49) 
aioe: 

Similar to the development of Iwagaki and Sakai for 
>, the coefficients A(x), B(x) and C(x) which have been 
deduced do not yield an exact solution of (30). Comparison 
of (39) and (49) shows that discrepancies exist in the third 
term coefficients of cos 2¢t and sin 2¢t (i.e., 6/5 versus 
1.0). The three terms which comprise the coefficients of 

-2 -5/2 8 


eos 29t and sin 2¢t involve x _, x and x ~. Use of 


= au and h = i*x enables the ratio of the third term 
to the first term to be evaluated as 
Rion: 
3rd term/lst term ~™ 2m (ah) (50) 
It is noted that the first and third term coefficients are 
in phase. For the case of i = 1/10, h = 25 cm and T = 3 
seconds, (50) indicates that the ratio of the third term 


to the first is equal to 0.089. Hence, the difference 


a) 





between the coefficients created by the use of (46), (47) 


and (48) is negligible and (40) may therefore be evaluated 


as 


gr t)\= cos2ct : sv SiIn{2Olx) 
«) aati eeaHen 


+ SIN 20t4 a (3f\' "<2 cos|2P(x) | 
-%%A_ <P sin 2 P(x] 


a 
a9 zie 


serra x ; (51) 


Equation (51) can be simplified to 


On = BWV *eeosfaleteze (ar Vem 


+%/2 + tan <a (ai) ee] E (52) 


In order to make this simplification, it is necessary to 


assume that 
t yo, 
l Uo a +1] =1.0, (53) 


The applicable range of solution for this perturbation 
scheme, which is discussed in Section V.C., restricts the 
computations to regions where 20 (EFS aoe cr othe 
lower limit of this relation, i.e. 20 (F)? =O. 2s 


substituted into the left side of (53), a numerical result 


of 1.166 is obtained. This error appears to be rather 
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significant. Utilizing a limit of 2¢ (2) 424.0, the left 
side of (53) is evaluated as 1.044. Consequently, the use 
of this ratio as a limit on the range of applicability of 
the solution may seem appropriate. Further discussion of 


this parameter is presented in Section V.C. 


B. SECOND ORDER EQUATIONS FOR 77 AND u 
Originally, it was assumed that 77 and u could be 
expressed as 


70 = 170 + ke 7 


(2) 5 


and 


(2) 


iu =o LO 40? Uae. = 


It is therefore necessary to combine the expressions for 
ny) and 7\) and for ys) and 62) to determine the 
final relations. It is noted that for (x,t), Iwagaki 
and Sakai neglected the second part of (22) which is inde- 
pendent of t. They felt justified in so doing since these 
terms effect only the stillwater depth. Their study was 
concerned with wave heights and profiles and consequently 
the stillwater depth was not required. This simplification 
is not valid for the present study since the phase speed, 
given by c = [er +h] 2, is effected by the stillwater 
depth. Furthermore, the ratio of the coefficients of the 


first term independent of t to the first term on the right 


DZ 





side of (22) is pa (£.)% i]. Utilizing the upper 
: 274 h 

limit for T CG) 2 of Gtr/i given by Iwagaki and Sakai, the 
importance of including these terms is apparent. (The upper 
icant of T (2)? is discussed in Section V.C.). 

When the first order equations for u and nr given by 
Iwagaki and Sakai, equations (20) and Cpe are expressed 
using the asymptotic approximations for the Bessel and 


Weber functions the following results are obtained, 


Seay: 
Wk t)= a | eo (ge) ae Tae 


+ 20 ( (ac - 1/4 (54) 


and 


U(x) a 0 ve. \2% CoS jot 
+20(ac)*— Ys 3. (55) 


Equations (22), (52), (54) and (55) can now be used to 


arrive at second order expressions for 77 and u, 


reatle al fe J EY 7a} | 


2 

GQ 
sof Dy cesta eet GG) *-= 
' Mae Ha ath 


+ = (% \e% cos e{ze(G) -% J 
= x * sinh be(z\* =H (56) 
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and 


u(x t)= «|= o( (4\4 "Gee ) eo sft 
meee WA § | 
eal $3 di 2 cos[2$ot+ 20 (gt VE} 
— eee “31 | (57) 


C. DETERMINATION OF THE PERTURBATION PARAMETER, 

The next question which must be addressed is how should 
the perturbation parameter, ~ , be defined. The first step 
in this determination, as outlined by Iwagaki and Sakai, is 
to formulate a region of applicability for the solution. 

In their study, they begin by considering a wave celerity 

of c = ne This is cause for some concern in that it 
differs from the accepted expression of c = [.g(77 +h)] * 

The problem, however, is to find a region of water where the 
theory can be applied; hence, use of an average velocity 

for the entire area seems appropriate. Prior to the arrival 
of the wave form and after its passage, the surface is at 
the stillwater depth and 77= 0. Therefore, 7( = 0 provides 
an average value for 77] over the region and c = (gh) 2 gives 
a representative average speed. In addition, over a large 
portion of the region 77 will be small compared to h. Ayyar's 


(1972) calculations indicate that even at the extreme point 
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of breaking, the phase speed at the crest will only vary 
from this value by a factor of (2)7, 
Using this value for c, Iwagaki and Sakai derive the 


relation 


| . 
h/AL = ca Ca (58) 


The authors conclude from this that an upper limit on h/L 
gives a lower limit for T(g/h)@. They then considered 
establishing an upper limit for T(e/h) 2. For this purpose 
they recalled that the approximations used for the Bessel 


and Neumann functions required that 





ee, 
26 (ar)* = Ie. (59) 


This in turn, implies that 


T(8h)* 4 ah Aa (60) 


Thus, (60) defines the upper limit of ieee Figure 5 
shows how (58) and (60) are combined to define the region 
of applicability for the solution. 

Iwagaki and Sakai select as the perturbation parameter 


the ratio hy/L > where h, is the largest depth in the 


i 


applicable range of the solution and Lo is the deep water 


> 
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Figure 5 
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wave length given by ce = oT /2ay . Use of this ratio may 
at first seem rather arbitrary. The motivation for this 
choice becomes more apparent when the derivation of the 
long wave equations by Freidrichs is recalled. His per- 
turbation parameter involved a representative length and 
depth quantity for the wave. Consequently, a similar ratio 
seems logical when considering the small value of &% . 

The ease with which this ratio may be evaluated pro- 
vides an additional incentive for its selection; examina- 
tion of the range of applicability of the solution is all 
that is required to obtain its value. An additional 
comment about this region is required. The upper limit 
to be placed upon h/L has not been previously discussed. 
In fact, two values for this parameter are shown in Figure 
5. For their numerical calculations, Iwagaki and Sakai 
use h/L €£1/20. The selection was arbitrary but conforms 
to general usage in wave theory. The theory was developed 
for shallow water where the pressure is hydrostatic; 
for this region, h/L £ 1/20 provides a reasonable limit. 


When this value is assumed, (58) requires that 


T (3/4 a (61) 


Applying Lj ,* qT Za to (61), 
h/t, =O @© la re (62) 
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Equation (62) will be adopted for the value of © throughout 
the majority of this presentation. 

Iwagaki and Sakai discuss the fact that use of (62) 
restricts the discussion to waves for which H/L, = 0.006. 
For greater Ho/Ly values, the theoretical energy flux curves 
predict that the deep water waves will break prior to 
arriving at the point where hy/L, =0-015/,-.1n order to 
extend the theory to situations where Ho/h exceeds 0.006, 
it is necessary to increase the upper limit placed upon h/L. 
Several cases for which h/L 41/15 are investigated in this 
study to determine the applicability of the theory to a 
domain of Ho/h values greater than 0.006. 

A final comment may now be made concerning the approxima- 
tion of (53). Use of (59) evaluates (60) as 1.16, not the 
assumed value 1.0. However, when these values are multiplied 
by the quantity ae as is required by (56) and (57), the 
difference becomes negligibly small. Hence the determina- 
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ae LO 





tion by Iwagaki and Sakai for the limit 2¢ ( Fi 


g 
seems appropriate. 
D. SECOND ORDER EXPRESSION FOR c 


A second order expression for the phase speed c can be 


found by substituting (56) into the phase speed relation, 


c= [en +n] %, giving 
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<= feta QQ ensfereetp 
ma | 


+ se % cosfa[etvae Gc) 
+ /> + tan oe Bo “eT 


ZS 
a ia *eosg2 [an (gi yeay 


Sa xP smfelen Gi)? 8 | 


GT 
iG 
a atx | = (63) 
The negative square root is utilized due to the fact that 


the wave is progressing in the negative x-direction. 


E. DETERMINATION OF a 

Prior to applying the kinematic breaking criterion to 
the preceding equations, the value of the term a in (56), 
(57) and (63) must be determined. This evaluation follows 
closely that which was outlined by Iwagaki and Sakai. As 
mentioned previously, these authors neglected those terms 
in a which were independent of t. Although these terms 
are included for the determination of c, they can be elim- 
inated for the purpose of establishing the value of a. The 
validity of this simplification stems from the fact that 


the value of the constant will be evaluated at the point 


Do 





where h 


=h 


ile the deepest depth in the applicable range of 


the solution. Examination of the terms independent of time 


in (56) shows that they decrease in absolute value for 


increasing x; the terms becoming negligibly small at the 


point where h = h,. Consequently, for the determination 


of a, 7 will be assumed to be expressed by 


7 (x it} ov [af | Ga} ha, sot 20(x “x4 | 


iE 
+ HK 2 |e ar x lcosjz[ot+2e(2e | 


: ae 


Substituting “ = hy/L, into (64) yields 


where 


r/h, = cere oo (26+ | (65) 


Ara Za i (1, ALY *(h, Jn) “(apn 
Pre ar ( hy.) (h, / h)lo/t J 


a= (ene) + afb. IE) po 
A. 


$= 1A + ten” SH (ony @ (b/A.) a) 
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1) 2. 
T£ AN and At ) are used to denote ao) and a(2) at 
h = h,, then from (66) 


"4 -3fy / 
Arar ha Ai (LAY a Hn, 


Agen! (n/L) (a/rn.) ) (67) 


In addition, at h = hy; panies (= Cees i(h/L,) 3) 
becomes negligibly small compared to t/2 and can therefore 
be neglected. 

The wave profile at h = h,, represented by ” 1? is 


now given by, 


2, 7h = A‘? ¢(e) (68) 


where 


F(@)= es oO ae SINCE, 


pa APA , Soe. 


At h= hy; the wave height is assumed to be twice the 


amplitude, 77 >? hence, 


H, /h,=2APE(&) (70) 


where @&c is evaluated from df/d& = 0. Iwagaki and Sakai 


determine Oc as 


Ol 





Cae ORC SIN ; Pit - ( 162)+2.) Vie G (71) 


The problem now is to find a value for H,/h,. Use of 
the theoretical curves for wave height change, which are 
based upon wave energy flux, provides a value for Hy /H, 
when h,/L, is known. H,/L,, the ratio of the deep water 
wave height to the deep water wave length, can be deter- 
mined for various wave conditions and then H, /h, is found 


through the identity 


H, /h, = CH, /H,)(Ho AAG /.,) we) 


Equations (67), (69), (71) and (72) can now be substituted 
into (70) to evaluate a. Simple computer techniques, as 
discussed in Appendix B, provide a determination of a for 


specific values of i, Ho/L, and period T. 
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VI. BREAKING CRITERION DERIVED 


A. APPLICATION OF KINEMATIC BREAKING CRITERION 

The derived horizontal wave induced velocity and wave 
phase speed are used to derive a breaking criterion. As 
stipulated by the kinematic breaking criterion, breaking 
will occur when the horizontal particle velocity equals 
the phase speed velocity. Theoretically, this condition 
may exist for several points in space and time. Of con- 
cern, however, is the specific case for which the horizontal 
distance from the beach at which u equals c is maximized. 
This will be the first position at which the approaching 
wave may break and hence all other cases are purely 
imaginary. The numerical solution technique employed is 
relatively simple. Subtracting the second order relation 
for u, (57), from that for the phase speed c, (63), equating 
the resulting expression to zero, and solving for x and t 
gives points in space and time at which the kinematic 
breaking criterion is satisfied. Examination of this 
solution set yields the maximum horizontal distance at 
which u equals c. The computer techniques employed are 


discussed in Appendix B. 
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B. RESULTS 

The numerical results obtained are summarized in Table 1. 
Several wave conditions are investigated in which the values 
for the beach slope, wave period, and deep water wave height 
to deep water wave length ey =) are varied. The specific 
selection of 0.119298 for the beach slope and 8.6 seconds 
for the period was made to conform with a future study 
(Hulstrand, 1976) in which experimental data will be used 
to verify these theoretical results. All other choices 
are strictly arbitrary. 

The first result of interest is that of the ratio 
ae Un» where 7, is the free surface elevation at the 
point of breaking. Table 1 shows that for each combination 
of beach slope and H/L, in the second order solution, the 
ratio is essentially constant (some small variations occur 
in the second significant figure). The value of the ratio 
is independent of the wave period. The consistency of 
the ratio suggests the use of this parameter as a breaking 
criteria for specific beach slopes and Ho/L, conditions. 

A second parameter listed in Table 1, that of h,/H,» 
where hy is the depth at breaking referenced to the still- 
water level, has often been utilized in the measurement of 
breaking waves. The theory under investigation yields 


consistent results for this ratio for specific beach slopes 
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and oe combinations. Mei (1966) derived a first order 


solution for this parameter of 


h/t, ~ (32)? (HA) a 


Examination of Figure 6 shows that Mei's calculations do 
not compare favorably with experimental data. Mei 
attributed this difference to the fact that his solution 
was confined to first order theory. Figure 6 shows that 
the theory presented here yields a much closer approxima- 
tion to the observed experimental data. It is noted that 
much of the observed data has been accumulated for values 
of H,/L, greater than those applicable to the current 
calculations. Plots for h/L & 1/20 and h/L <1/15 have 
been extended through the theoretically derived points for 
purposes of comparison. Within the concurrent regions of 
applicability, the theoretically derived values for the 
various limiting values of h/L are similar. This suggests 
that the extension of these graphs into the domain of 
greater Ee ee values provides at least an indication of the 
breaking criteria which would be derived by utilizing h/L 
ratios applicable for these regions. 

The theoretically derived values of hy /H, decrease 
with increasing values of Ey eee This conforms with the 


trend displayed by the field measurements. In contrast, 
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Mei's results show an increase in the ratio with increasing 
H,/L- 

Of concern is the fact the current theory predicts 
breaking depths less than those observed in the experimental 
data. Table 1 and Figure 6 include data for the first 
order solution resulting from the theory presented. Figure 
6 shows that the extension to the second order solution pro- 
duces results which compare more favorably to the experi- 
mental results. This suggests that higher order solutions 
would yield improved results. An additional source of 
error present in this study results from the fact that the 
vertical water particle accelerations have been neglected. 

A final comment can be made concerning Mei's calcula- 
tions. The first order solution which is derived from the 
procedure utilized in this study yields considerably 
improved results over those of Mei. Therefore Mei's 
descrepancies cannot be attributed entirely to the restric- 
tion to the lowest order solution. The source of error may 
possibly be due to the use of the geometrical breaking 
eriterion. 

Figures 7 and 8 depict representative wave profiles at 
breaking for the current theory. These indicate that 
breaking occurs prior to the attainment of a vertical 


surface slope. The figures also show that the predicted 
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breaking occurs at the wave crest. Figure 8 shows an 
unrealistically deep trough in front of the crest. The 
use of the Bessel function approximation has been shown 


Saiic 


to be valid for 20(5 2 


1.0 . Substitution into this 
relation of the values of the terms used to construct 
Figure 8 requires x221.89 feet. It is seen, chereteee 
that the excessively deep trough is predicted in a region 
where the Bessel function approximation is not valid. The 
predicted breaking for the specific instance shown in the 
figure is at x = 33.73 feet, which is in the applicable 
region for the approximation. For each case investigated, 
the predicted breaking point occurred well within the 
region where the approximation is accurate; hence the 
derived breaking criteria is deemed valid. 

Ayyar (1970), utilizing the kinematic breaking criterion, 
produced the additional shallow water breaking criterion of 
Yp/h,=2.0. Table 1 summarizes the values obtained in the 
present study for this ratio, all of which are considerably 
larger than 2.0. Ayyar's ratio does not account for the 
variations associated with Hy/L, suggested by the observed 
h,/H, data. In addition, his theory is limited to waves 
which have the geometry of a plunging breaker. Ayyar also 
assumes that breaking occurs at the wave front, which may 


not in fact take place. The limitations placed upon Ayyar's 





formulation and the random values of the ratio displayed 
in Table 1 prevents this from being considered a valid 
indacation of breaking. 

The most significant disadvantage associated with the 
breaking criteria derived from the theory presented in 
this study is that the ratios require re-calculation for 
changes in beach slope and deep water wave height to wave 
length ratio. This objection is similar to that which 
was raised against Stoker's use of the method of charac- 
teristics. The complex dependence of the values of a, 

u and c upon both the beach slope and H/Uy: however, 
leaves the investigator with little choice but to resort 
to a numerical solution. In defense of the approach, the 
solution technique applied to specific situations is 
relatively simple and requires minimal computations once 
the beach slope and Ho/k, are known. pellecced data points 
may be used to construct graphs which approximate the 

i, H ratio for each specific beach slope. This can be 
used to provide an indication of breaking for varying 


ele values. 
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VII. CONCLUSIONS 


Two shallow water breaking criteria have been formu- 
lated through the application of the kinematic breaking 
criterion to a second order solution of the long wave 
equations. These are the ratios of 7,/L, and h,/H,. 

Both of these criteria are dependent only upon the beach 
slope and the ratio of the deep water wave height to the 
deep water wave length. Each ratio requires re-calculation 
as these two parameters vary. Comparison with previous 
theory indicates that the theoretically derived values 

for hp/H, offer significantly improved approximations to 
the assembled experimental data. The increased accuracy 

is partially attributed to the use of a higher order 
solution to the long wave equations. First order solutions 
obtained indicate that the solution technique applied 
offers improvement over previous theory. The predicted 
breaking depths are somewhat less than those observed in 
field measurements. It is believed that the extension of 
the theory to higher order solutions of the long wave 
equations would reduce the error associated with the pre- 
dicted breaking depth. Inclusion of vertical water particle 


accelerations would also increase the accuracy of the solution. 





APPENDIX A 


DERIVATION OF 7 (2) 


Iwagaki and Sakai have obtained a second order solution 
for the free surface, (x,t). The relation derived com- 
pared favorably with experimental data. Several approxi- 
mations made in the solution, however, warrant discussion. 

Combining equations (18) and (19) so that 62) 


eliminated, gives 
ee. - gl np ix +7, | Tis 
) 
- [> a] . teu Pix + Lut re 
Iwagaki and Sakai substituted their first order equations 
(20) and (21) into the right side of (A). They then 


offered approximations for the Bessel and Neumann (Weber) 


functions of 


J-(w)™ (26) *%e os Cw ay) aa 1 /4 ) (A-2) 


and 


Ne Cw) (hin)! "910 (w- rV/2- 1/4 ). (A-3) 


The authors provide evidence which shows that for lw] >1.0, 


these asymptotic expansions are accurate. Defining 





§ (x) = (2 ox (x/gi)2- 71/4), use of (A-2) and (A-3) allows 


the right side of (A-1) to be expressed as 


pens { 308 2)" “cos| 2 P(x) | 
a % at os Bia A oul 
+ 2a iy x Zea s t3 PIE 


tsmdot § 38 (Bc % sm [2G (0) 
ef [zo ish a ogee (x 5 / 
oe ee x 4 SIN on é 


+f oe (BY* x % cos EH] 
va es Asin) 2 I(x) ] 


~| 


- a+ (at) ye ee s| 2 GX] : 


(A-4) 


The solution of > was assumed ta be 
AX) =A(x) cosfal+ B(x) sin Got + CC). 
This expression was substituted into the left side of 
equation (A-1) and the result compared ‘to (A-4). From this 
comparison, the coefficients A(x), B(x) and C(x) are deter- 


mined as 


ae te (a\/* aks | 2 P(x], (A-5) 








B(x) 3 a x cos |2pP (x) | 
as x se sin |] 2P(x) | (a6) 
C (x)= 4 z= ( a “72 cos[Z2P(x) | 


Ss ao = SiN \29 (x) }. (A-7) 


Use of (A-5), (A-6) and (A-7) gives the left side of (A-1) 


as 


ae ce Bae (BNA -¥2 os [2H(x) | 
“fi piaa: “3m a eel29 Od 
+20 a4 (aii cosle Gul 


+ SIN zat § ae x ae aera O(x)] 
+ Sha gfrx* es Pla) ] 
ee og (404: sin 2 Pla] § 


mee 

+L acfy CBR t0s)ZP (x) { 

ee ae 3 /ar x 2 LEA | 

Ts (gt) x eae PC x) } 

rl 29, x Sein) POOL T. (A-8) 
Examination of (A-4) and (A-8) shows that use of the 
expressions for the coefficients A(x), B(x) and C(x) does 
not yield an exact solution. Differences occur in the 
third term constants in each of the coefficients of 


cos 2¢0t and sin 2¢0t. In addition, the segment of (A-8) 


independent of t contains a fourth term not present in (A-4). 


ray 





Instinctively, these differences would seem to limit 
the accuracy of the solution. Iwagaki and Sakai substan- 
tiated their solution, however, by considering the relative 
Significance of the terms. It was noted that the terms of 
the coefficients of sin 27t and cos 2¢t involved the 


5) 2 5/2 


values x ~'°, x * and x ~/“, while the fourth term inde- 


pendent of t in (36) contained x7>, Using og = 27/T and 
comparing the second, third and fourth terms to the first, 


the following ratios were found, 


Padsberm/ist Term ~ be { (aT zr) Ks L) as 
3rd Term/lst Term ~ [ 3 (qrZn)/( Lye Lo 


and 


4th Term/lst Term ~ Ei (qT ew Wary te 


The interpretation of these ratios is that the successively 


higher terms become relatively smaller in proportion to i. 
i 
10 ° 


h = 20 cm and T = 3 sec. Examination of the first and third 


Iwagaki and Sakai considered the specific case of i = 


terms show that they are in phase. For these particular 
conditions, the ratio between the third and first terms is 
less than =: The conclusion is that the difference 


between the constants for the third terms is negligible. 
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Similarly, the fourth term is compared to the second 


which is in phase with it, and seen to be negligible. The 


use of (A-5), (A-6) and (A-7) can therefore be used to 


e 2 2 
provide an accurate expression for % ) as 


7x 6)= cgs2at 3- oma Ces o(x) | 
vo ae (2\45 2 oslL2 Pin] § 


+ SIN 20t§ = a * <os]2O(x) | 
~ "0 $e ae (ali) 7% sin T2 ly) | § 
+] ak (8, yes -42.cos 129 ] 
~ Ke og sin [290] J. (A-9) 


The final solution for fh?) (x,t) offered by Iwagaki 


and Sakai was 
7O°Xx Et= oe oe Ne cos} 2[ot r20(aeyM/4] 
ead (x S aio™ = 
fe = Ce (2 x 2 <2 sfoceate ye vA) 
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APPENDIX B 


COMPUTER PROGRAMS 


The results which are summarized in Table 1 were calcu- 
lated through the use of a series of simple computer pro- 
grams. They are essentially a series of do-loops designed 
to perform a number of iterations over an interval. The 
results are examined to determine the desired solutions. 

The first two programs used are concerned with estab- 
lishing a value for the quantity, a, found in equations 
(56), (57) and (63). As was discussed in Section V.E., 
this evaluation can be made through the use of (67), (69), 
(71), (72) and (70). The first step in the process is to 


employ (72), 


Hy /hy i OSG Liebe) CsCl) 7 One ees (72) 


Eo, Lind Hy /h,. In this expression, hy/L, is known for the 
particular case under consideration. Equation (62) gives 
this term as h,/L, = 0.0157 for the limiting condition of 
h/L £4 1/20. Entry into the wave energy flux curves of 
hyperbolic wave theory with the value of h,/L = 0.0157 
gives a value of H, /H, = 1.28. The remaining quantity in 


the right side of (72), H/Lo, is a known value for the 





specific wave conditions. It can be determined from deep 
water observations of the wave. Therefore, once this 
value is specified for a certain set of wave conditions, 
(72) yields a value of H,/h,. For instance, when Ho/L, 
is equal to 0.001, H,/h, amino) C001) ) 20157 = 0, 0815. 
Equation (69) provides a relation for the quantity b 


2 
as b = at yaQ), Use of (67), which defines al?) and 


yb) 


1? allows this equation to be evaluated as 


b= 74 2M (nih) *ah,, ; (an) 


Z 
Bee ol = 0.0157 and L = 2! (B-1) is equivalent to 
E 0 Zot 


oe nits ae Gra gt @0IF7)* (B-2) 


The next value which must be determined is that of ECC) 
in (70) which is defined in (69) as 
F(G@.) = cos Oe —bSIN 20. - 


Use of the identity 


Sir 20= _2ten GO 


ee ano 


reduces (69) to 


F(G)= cos Og — aS as ‘ (B-3) 





Equation (71) gives Oc as 


es ace sm §(de- Gig 4)*)/28. 


Defining x as, 


ae ie (B-4) 


and using the identities which exist when «@ = arc sin x of 


Cos eC = Jol Sac 

and 
| Zz x (B-5) 
ee 


results in G3) being given by 


F(@,)= | [ino eee = 


a a B-6 
i+ (=) (B-6) 
(1). 
A, is determined from (67) as 


ee ee ee, Nc” - 
5 oie Ch, Co0'57)aT~ (B-7) 


Equations (B-6) and (B-7) can now be combined to express 
Hy /hy in terms of the quantity a. The requirement then is 
that of determining the value of a for which the right side 
of (70) equals the value of H,/h, as given by (72). 

The first computer program used defines a function F 


as the difference of the right hand value of (70), obtained 
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through (B-6) and (B-7), and the left side evaluated by 
(72). The program is designed to plot F as a function of 
a. The value of a for which F equals zero is the required 
solution of (70). Examination of the graph produced by 
the first program can thus give a first approximation to 
the quantity a. 

The second computer program employed is simply an 
iteration routine used to refine the value of a. It begins 
with a first estimate of a determined from the graph of F 
versus a. The value of a is then incremented in steps of 
0.001 and the corresponding values of F calculated. The 
program is designed to determine the value of F closest to 
zero and to print the value of a for this case. This is 
the value of a which satisfied (70). 

The next series of programs are designed to determine 
the maximum value of x at which the kinematic breaking 
criterion is satisfied and the time at which this occurs. 
This can be accomplished by finding a solution set of (x,t) 
combinations for which the difference between the right sides 
of equations (63) and (57) is equal to zero. The technique 
used is similar to that used in establishing the value of 
a. The presence of two dependent variables, however, 


makes the procedure somewhat more complex. 
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The first program in this determination utilizes the 
method employed for the plot of F as a function of a. A 
function is defined which equals the difference between 
(63) and (57), (c-u). This is plotted as a function of x 
for a series of times t. Examination of the graphs provides 
-a first guess at the maximum x for which the kinematic 
breaking condition is satisfied and the time at which it 
occurs. As was the case for determining a, a second pro- 
gram is now utilized to obtain a more refined solution. 
Here, time is varied over one second in steps of 0.05 
seconds and x is incremented in intervals of 0.01 feet. 
The value of the function defined as (c-u) is printed for 
these specific (x,t) combinations. These results are 
examined to determine the point at which the kinematic 
criterion is first satisfied. 

These first programs have thus determined the maximum 
x at which the kinematic breaking criterion is satisfied 
and the time at which it occurs. The remaining programs 
utilize this result to produce the breaking criteria 
summarized in Table 1. The fifth program in the complete 
series simply calculates the value of 77 (the free surface) 
and the depth of the bottom for the point specified for 
the solution of (x,t). This is accomplished by using 


(56) to define 77 and h = i*x,. The ratios of p/h, and 
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and h,/L, are then calculated. A similar program is also 
provided which graphs a as a function of x. 
Pages 86-91 contain flow charts for the programs 


employed in these determinations. 
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PLOT 
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STOP 


PROGRAM 1: APPROXIMATION OF a 


Figure 9 
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Figure ll 
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